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ABSTRACT 

During the late stages of their evolution, Sun-like stars bring the products of nuclear 
burning to the surface. Most of the carbon in the Universe is believed to originate 
from stars with masses up to a few solar masses. Although there is a chemical di¬ 
chotomy between oxygen-rich and carbon-rich evolved stars, the dredge-up itself has 
never been directly observed. In the last three decades, however, a few stars have 
been shown to display both carbon- and oxygen-rich material in their circumstellar 
envelopes. Two models have been proposed to explain this dual chemistry: one pos¬ 
tulates that a recent dredge-up of carbon produced by nucleosynthesis inside the star 
during the Asymptotic Giant Branch changed the surface chemistry of the star. The 
other model postulates that oxygen-rich material exists in stable keplerian rotation 
around the central star. The two models make contradictory, testable, predictions on 
the location of the oxygen-rich material, either located further from the star than the 
carbon-rich gas, or very close to the star in a stable disk. Using the Faint Object In- 
fraRed CAmera (FORCAST) instrument on board the Stratospheric Observatory for 
Infrared Astronomy (SOFIA) Telescope, we obtained images of the carbon-rich plan¬ 
etary nebula (PN) BD -1-30° 3639 which trace both carbon-rich polycyclic aromatic 
hydrocarbons (PAHs) and oxygen-rich silicate dust. With the superior spectral cover¬ 
age of SOFIA, and using a 3D photoionisation and dust radiative transfer model we 
prove that the 0-rich material is distributed in a shell in the outer parts of the nebula, 
while the C-rich material is located in the inner parts of the nebula. These observations 
combined with the model, suggest a recent change in stellar surface composition for 
the double chemistry in this object. This is evidence for dredge-up occurring r\j 10^ yr 
ago. 
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1 INTRODUCTION 

Planetary Nebulae (PNe) are the final evolutionary phase of 
low- and intermediate-mass stars. The nebulae form out of 
the mass lost by the star on the asymptotic giant branch 
(AGB), which may exceed 50% of the stellar mass. The 
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ejecta mainly consist of gas, initially molecular and atomic 
but becoming ionised by the remnant white dwarf star. Some 
solid particles (‘dust’) condense in the ejecta. The ejecta 
quickly disperses and merges with the surrounding inter¬ 
stellar medium. This recycled gas fuels most of the star for¬ 
mation in late-type galaxies ( Leitner Kravtsov||20II ). 


The chemical evolution of the Universe is driven by 
products of nucleosynthesis included in stellar ejecta. Low- 
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mass stars are the primary source of new dust in the Milky 


Way (Matsuura et ah 2009), and the dominant producers 


of carbon and nitrogen ( Kobayashi et aL]|2011 Henning & 

Salama|1998 ) 


Low- and intermediate-mass stars produce PNe with a 
distinct molecular and dust composition. In the molecular 
zone, the highly stable but volatile carbon monoxide (CO) 
molecule locks away the less abundant element, leaving the 
remaining free oxygen (O) or carbon (C) to drive the chem¬ 
istry and dust formation. If C/0<1 all the carbon is trapped 
in CO, and the chemistry forms oxides and silicate dust, with 
spectral signatures at 9.8, 18.0, 23.5, 27.5, and 33.8pm. If 
C/0>1 then carbon dominates the chemistry and the main 
observed dust features are amorphous carbon (no spectral 
feature), SiC, complex hydrocarbons, including polycyclic 
aromatic hydrocarbons (PAHs) (features centred at 6 . 2 , 7.7, 
8 . 6 , and 11.3pm), and fullerenes ( |Cami et al.||201Q ). Such a 
dichotomy between O- and C-rich PNe is observed. Stars 
begin their life with C/0<1 and are thus oxygen-rich. In 
the AGB phase, the dredge-up of newly synthesised carbon 
during a phase of helium flashes can raise the C/0 ratio 
above unity to form a carbon star. The details of this process 
are still very uncertain (Karakas|2011| Karakas & Lattanzio 
20T4|. 


A small fraction of PNe show both O- and C-rich 
features in their dust spectra, and are therefore classi- 


fied as dual chemistry objects ([Waters et al.| 1998 

11 ^ 

stra et al. |1991| [Waters et al.[[1998| [Cohen et al. 

[1999 

2002 [ [Gutenkunst et al.[[2008 

Perea-Calderon et al. 

2009 

Guzman-Ramirez et al.[[2011| 

2014). Such objects can pro- 


vide indirect evidence for dredge-up and allow the study of 
the environment in which it takes place. Two models have 
been proposed to explain this dual chemistry. One model 
requires oxygen-rich material to be present in an old and 
stable disk around the carbon-rich central star, and is not 


related to recent mass loss (Zijlstra et al. 1991). There is 
evidence for the existence of old, stable disks around post- 


AGB stars (de Ruyter et al. 2006). The other postulates 


that a recent dredge-up of carbon produced by nucleosyn¬ 


thesis changed the surface chemistry of the star (Waters et 


al. 1998). Finding objects in the process of changing from 


0-rich to C-rich chemistry has remained elusive (Zijlstra et 
al.|2004| |. 


BD +30° 3639 (PNG 064.7+05.0) is one of the few plan¬ 
etary nebulae to host a Wolf-Rayet [WC] central stars. These 
objects show emission-line spectra similar to those of Popu¬ 
lation I Wolf-Rayet stars but have lower masses, expected for 


intermediate-mass stars in the post-AGB evolution (Leahy 
et al.|[20QQ ). Their spectra show strongly enhanced carbon 


and helium but have little or no hydrogen in the atmo¬ 
sphere. The central star of BD +30° 3639 (HD 184738) is a 
hydrogen-deficient, carbon-rich ([WC9]) star ( |Mendez|1991 ) 
with an effective temperature of 55,000 K, a lumi nosity of 
4.25x10^ L 0 , and a post-ejection mass of 0.6M© ( Li et al. 
|2002| . BD +30° 3639 is active in all spectral regions. While 
most of its energy is emitted in the mid-infrared (peaking at 
~30pm), it is also a strong radio source (Zijlstra et al. IToMl i. 
It also has a molecular envelope with a mass of 0.016 Mq es¬ 
timated from millimetre rotational lines of CO (Bachiller et 
al.|1991| ). 


2 OBSERVATIONS 

Using the Faint Object InfraRed CAmera (FORCAST) 
for the Stratospheric Observatory for Infrared Astronomy 
(SOFIA) Telescope ( Adams et al.| 20 l 0 ) we obtained images 
of BD +30° 3639. The observations cover the wavelength 
ranges of both the PAHs and crystalline silicate features. 
SOFIA uniquely allows for the imaging of the long wave¬ 
length emission from crystalline silicate features at high an¬ 
gular resolution. 

The FORCAST instrument has a short wavelength 
camera (SWC) that operates from 5-25pm and a long wave¬ 
length camera (LWC) that operates from 25-40pm, with 
several filters available in both cameras. FORCAST sam¬ 
ples at 0.75'Ypixel, giving a 3.2'x 3.2' instantaneous field- 
of-view. Five different filters were used: FOR_F064 (cen¬ 
tred at 6.4pm with a AA=0.14pm), FOR_F077 (centred at 
7.7pm with a AA=0.47pm), FOR_Flll (centred at 11.1pm 
with a AA=0.95pm), FOR_F113 (centred at 11.3pm with 
a AA=0.24pm), and FOR_F336 (centred at 33.6pm with 
a AA= 1.9pm). We used the symmetric nod-match-chop 
(NMC) imaging mode, the chop is symmetric about the op¬ 
tical axis of the telescope with one of the two chop positions 
centred on the target. The nod throw is oriented 180° from 
the chop, such that when the telescope nods, the source is 
located in the opposite chop position. 


3 RESULTS 

In Figure 1 we present the images of BD +30° 3639 obtained 
using the filters centred at 6.4, 7.7, 11.1, 11.3, and 33.6pm. 
Each image has an over-plotted blue circle of 5" in radius. 
To study the spatial extension of the different dust features 
in the PN, we performed a radial cut along the major axis of 
the nebula (PA=135°) for all the SOFIA images (the white 
line in the top left panel represents the PA of this cut). Fig¬ 
ure 2 shows these cuts, the top plot shows the cut made 
using the PAH filters (centred at 6.4pm, 7.7pm, 11.1pm, 
11.3 pm). The extent of the shells in the PAHs filters is very 
similar. The cuts of the 7.7 and the 11.3pm images show a 
double peaked emission, as expected for a shell surrounding 
a hollow interior. The measured full width at half maxi¬ 
mum (FWHM) is 8 " at these wavelengths. The middle im¬ 
age shows the PAHs filters and the silicates filter (centred 
at 33.6pm). The emission at 33.6pm is 1.5 times more ex¬ 
tended than the C-rich material, with a FWHM of 12". The 
bottom plot shows the same radial cuts done in all the filters 
with the normalised flux. 

Figure 3 shows the Infrared Space Observatory 
(ISO) Short Wavelength Spectrometer (SWS) spectrum of 
BD +30° 3639 (reproduced from Sloan et al.|2003| ), with the 
bandpass of the SOFIA filters overlaid by a dashed line of 
different colours depending on the filter. Dark-blue repre¬ 
sents the bandpass for the FOR_F064 filter, the FOR_F077 
filter bandpass is in purple, brown for the FOR_Flll fil¬ 
ter, green for the FOR_F113 filter, and light-blue for the 
FOR_F336 filter. Using the ISO spectrum we can clearly see 
the PAHs and the silicates features covered by each filter. 
We can see the main PAHs (6.25pm, 7.7pm, and 11.3pm) 
as well as the crystalline silicates (33.6pm) features. 

Previous observations of BD +30° 3639 show a shell-like 
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structure with a radi us of 5" x A" in the visible (Li et ah 2002) 
and in radio ( Bryce et al.|l997 ), while a much larger 

extension of the faintest part of the nebula is detected in 
the near-infrared ( Phillips Ramos-Larios||2007 ). Previous 
mid-infrared imaging observations of BD +30° 3639 with a 
spatial resolution of 1.5'' found that the 8.6 and 11.3pm 
(PAHs) bands are slightly more extended than the contin¬ 
uum emission ( Bentley et al.|19^ ). Higher spatial resolution 
imaging and spectroscopic observations of BD +30° 3639, 
taken with the Subaru telescope, confirmed that the PAH 
emission bands are more extended than the dust continuum 
emission, but show a similar extent to the [Ne H] 12.8pm 
emission ( Matsumoto et a'r]|2008 ). These observations also 
show evidence for a 10 pm absorption feature, attributed to 
silicates. 

Until SOFIA, it was not possible to map the silicate 
emission features in this object, because the features emit at 
longer infrared wavelengths where no instrument was able to 
offer the necessary combination of wavelength coverage, sen¬ 
sitivity, and spatial resolution. Using the expansion velocity 
of the silicate shell of 10 km/s, and a distance of 1.2 kpc 
( Li et al.|2002 ), the expansion rate of the shell is 1.76+0.25 
masyr“^ in radius. Using the images we compared the size 
of the C-rich to the 0-rich shell. For the C-rich shell we av¬ 
eraged all the images of the PAHs filters (6.4 to 11.3pm), 
for the 0 -rich shell we used the image of the silicates filter 
(33.6pm). We measure the size (radius) of the C-rich and 
the 0-rich shells to be 5+0.75" and 7.5+0.75", respectively. 
These correspond to dynamical ages of 2,800+580 yr for the 
C-rich shell and 4,300+740 yr for its 0-rich counterpart. 
This means that the transition of the star from O- to C-rich 
took place within a window of 1,500+940 yr. These values 
are entirely consistent with the post-AGB age of the central 
star determined by comparison of its luminosity and tem¬ 


perature to evolutionary tracks (Vassiliadis & Wood 1993 
BlQcker||2001 ). 


4 ANALYSIS ^ DISCUSSION 

In Figure 1 and 2, we can see the different extensions of the 
two major components presented here, PAHs and silicates. 
It is important to note that no continuum filter observations 
were made, as such, the dust continuum contribution is not 
subtracted from any of the presented observations. There¬ 
fore, one must consider the level of continuum contamination 
in each filter before drawing conclusions as to the relative 
distribution of PAHs and silicates. 

In the case of PAHs, from Figure 3, we can see that the 
FOR_F064 filter has a negligible contamination by dust con¬ 
tinuum and thus traces the PAHs emission. From Figure 2 
we can see that the emission in this hlter has the same ex¬ 
tension as the other PAHs filters, as such, we can assume 
that the measurements made in the PAHs filters are tracing 
the PAH emission, rather than just extended dust contin¬ 
uum. 

Unfortunately, the observation made in the filter encompass¬ 
ing the silicates emission, the FOR_F336 filter, includes a 
significant contribution of dust continuum, as shown in Fig¬ 
ure 3. Based on the images alone, we are unable to determine 
the spatial distribution of the silicates with respect to the 
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Figure 1. SOFIA images of BD+30° 3639 using the 6.4pm, 
7.7pm, 11.1pm, 11.3pm, and 33.6pm filters. To compare the size 
of each image, we have overplotted a blue circle of 5" radius in 
all of them. The white line in the top left panel represents the PA 
of the cut made to all the images. 


dust continuum in the nebula. In order to resolve this degen¬ 
eracy, a radiative transfer model of the nebula was created. 
We constructed simple 3D photoionisation and dust ra¬ 


diative transfer models using the MOCASSIN code (Ercolano 


et al. 


2003 


2005 


2008). The model consists of neutral 


PAHs and graphitic carbon (optical constants from 


Draine 2001), and amorphous silicates (optical constants 
from Laor Draine||1993 ) and crystalline fosterite (optical 
constants from Jager et al.||1998 ). In both the dust density 
is proportional to r“^ and the grain size distribution is 
standard Mathis, Rumpl, & Nordsieck (1977 [MRN]). We 
adopted a distance of 1.2 kpc, and a central star luminosity 
and temperature of 4.25xlO^Lo and 55,000 K. We varied 
the shell radii, the position of the silicates and carbon, 
and the mass of dust in each shell to obtain a good fit 
to the data. For the model to fit the ISO spectrum and 
the SOFIA data points, the carbon material has to be in 
the inner parts, while the silicates need to be in the outer 
parts of the nebula. If silicate dust is present in the inner 
warmer regions of the shell, emission features at 10 and 
18pm are predicted, opposite to what it is observed. The 
total mass is constrained to within 10% by the models. In 
the best fitting model (Figure 4), the carbon shell extends 
from 1.5 to 5.5" on the plane of the sky and contains 
6 . 6 xlO~^M 0 of dust, while the silicate shell extends from 
5.5 to 7.5" and has a mass of 8.1x10“^ M©. The carbon 
dust mass loss rate, assuming the switch happened 2,800 yr 
ago, would be 2.4x 10“® M© yr“^, so a total mass loss rate 
of 7.0x 10“® M© yr“^ if the gas-to-dust ratio is 300. The 
silicate dust mass loss rate would be 5.4x10“® M©yr“^ for 
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Figure 2. Radial cuts along the major axis of BD +30° 3639, 
along a PA of 135° for all the SOFIA images: Top image shows 
the PAHs filters only (6.4|xm, 7.7|xm, ll.l|xm, 11.3|xm). The mid¬ 
dle image shows the PAHs filters and the silicates filter (33.6|J.m), 
while the bottom image shows all the filters with the flux nor¬ 
malised. 


FOR_Flll 

FOR_F064 FOR_F077 FOR_Fll3 FOR_F336 



Figure 3. ISO spectrum of BD +30° 3639. The dashed lines are 
used to show the position and width of the SOFIA filters used on 
this paper. Dark-blue represents the bandpass for the FOR_F064 
filter, the FOR_F077 filter bandpass is in purple, brown for the 
FOR_Flll filter, green for the FOR_F113 filter, and light-blue for 
the FOR_F336 filter. 



a total mass loss rate of 1.6 x 10 ^ M© yr 


The observed C/0 ratio in the ionized region is 1.59 
( Bernard-Salas et al.|2003| ). Assuming instant, homogeneous 
mixing, and a solar O abundance, an initial C/0 ratio of 0.8, 
and an envelope mass of O.OI M©, the numbers suggest that 
the thermal pulse dredged-up 4 x 10“^ M© of carbon, this is 
almost the same value we obtain. However, the current wind 
from the central star has C/0 = 30 zb 15 ( Yu et alT||2QQ9 ), 
so that even if the mass-loss rates are low, the current wind 
can still add significantly to the carbon budget. The range of 
C/0 and the somewhat low carbon dust mass may suggest 
that the carbon dredge-up was not instantaneous but that 
the C/0 increased while the mass loss continued. 

Independently, a stellar evolutionary model was con¬ 
structed with an initial stellar mass of 1.5 M© and solar 
metallicity, fairly typical values for a Galactic Disk PN. This 


Figure 4. 3D model using MOCASSIN. The blue dotted line is the 
total SED given by the model, in cyan we show the contribution 
of the O-rich dust and in pink the C-rich dust. The green crosses 
are the SOFIA fluxes (assuming a 20% uncertainty), and the red 
line is the ISO spectrum over-plotted. 


model (Figure 5) derives a peak mass-loss rate at the point 
where the model star becomes C-rich of ~ IxIO ^ M© yr ^ 


and a final central star remnant mass of 0.62 M© (Karakas 


2014), closely matching those values determined from the 


observations and photoionisation modelling. The model 
shows that the typical interpulse period near the end of the 
star’s life is ~IxI0^ years, while the total duration of a 
thermal pulse, a third dredge-up, and relaxation back to H- 
shell burning is of the order of ~I000 years or more. The 
numbers fit the duration of the change from oxygen-rich 
to carbon-rich mass loss observed with SOFIA, and con- 
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Figure 5. Evolution of the mass-loss rate with time, as estimated 
from our stellar evolution model. We used an initial mass of the 
star of 1.5Mq, and a solar metallicity (Z = 0.014). The peak 
mass-loss rate is ~lx 10“^ Mq yr“^ and the final mass of the 
star is 0.62Mq, which closely matches the value calculated from 
the observations. 


firm that this is the likely explanation. At the last thermal 
pulse, 4.0 X 10“^ Mq is dredged-up to the envelope. About 
1 X 10“^ Mq of this is This is higher than derived from 
the observations. Incomplete mixing (Hajduk et al. 2005) 
can therefore not be ruled out. 

The observations and modelling suggest that the ther¬ 
mal pulse took place while the envelope mass was already 
very low. This is consistent with an AGB final (or fatal) 
thermal pulse (AFTP) occurring immediately before the star 
moves off the AGB. An AFTP can lead to both a consider¬ 
able enrichment with carbon and oxygen and to the dilution 
of hydrogen, and also explains the current hydrogen-poor 
nature of the current central star ( Blocker [200 1 ). 

The short duration of this phase compared to the long 
interpulse period suggests that such observations of carbon 
dredge-up will remain uncommon. BD+30 3639 provides a 
rare glimpse into a phase of stellar evolution which is crucial 
to the origin of carbon in the Universe. 
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